AD-A051  863  INOIANA  UNIV  BLOOMINGTON  DEPT  OF  MATHEMATICS  F/6  12/1 


ASYMPTOTIC  MULTINORMALITY  AND  REMAINDER  TERMS  OF  LINEAR  RANK  VE~ETC(U) 
FEB  78  M L PURI*  N S RAJARAM  AF0SR-76-2927 

UNCLASSIFIED  AFOSR-TR-78-0434  NL 


AFOSR-TR-  7 8-  0434 


I j 


(l.D 


c2 


Asymptotic  Multinormality  and  Remainder  Terms  D D C 
of  Linear  Rank  Vectors  Under  Alternatives  t — vr-i/ — v. — inn  nr 

[n]E@!ZDn  f7(7 

bv  MAR  28  1973 


Madan  L.  Puri  and  Navaratna  S.  Rajaram 
Indiana  University  and  Kent  State  University 


0.  Summary 


Asymptotic  multinormality  of  linear  rank  statistics 
based  on  independent  vector  valued  random  variables 
is  obtained.  Under  suitable  assumptions,  weak  esti- 
mates for  the  remainder  terms  for  convergence  to  normality 
are  also  obtained.  Results  on  asymptotic  normality  are  re- 
lated to  Ha'jek  (1968)  and  Puri  and  Sen  (1969).  Re- 
sults on  the  remainder  terms  are  related  to  those  of 
Jureckova''  and  Puri  (1975),  Bergstrom  and  Puri  (1977),  and 
Huskov^  (1977)  . 


' |jlJ  1.  Preliminaries 


Let  = (X^]...,  X^?))  , i=l,...,N  be  a sequence  of 
independent  p-variate  (p  a 1)  random  vectors  having  con- 
tinuous cumulative  distribution  functions  (cdf)  F .(x), 

Ni  ~ 

5 = , . • . ) , i=l,  ...,N  respectively.  Consider 

now  the  random  matrix  XN  corresponding  to  ... # i,e' 

= ((XNi}))  : &N1  “ <xiii> XN?))' 

1- if  • « • f N| 
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ABSTRACT  (Continue  on  reveree  aide  If  neceaaary  end  Identify  by  block  number) 


Asymptotic  multinormality  of  linear  rank  statistics  based  on  inde- 
pendent vector  valued  random  variables  is  obtained,  under  suitabli 
assumptions,  weak  estimates  for  the  remainder  terms  for  convergent 
to  normality  are  -related  to  Hajek  (Ann.  Math.  Stat.  (1968V  3 2 5-346). 

Results  on  the  remainder  terms  are  related  to  those  of  Jureckova 
and  Puri  (Ann.  Prato.  (3975),  526-  533),  and  Bergstrom  and  Riri  (Ann.  Stat. 

(1977),  671-680)  . Ctt  ff  


DD  1473  EDITION  OF  1 NOV  «S  IS  OBSOLETE 


E5raS«= 


LAM!  FIC  ATION  OF  THIS  FACE 


-2- 


(1.2) 


(1.3) 


(1.4) 


(1.5) 


and  observe  that  each  row  of  is  composed  of  N inde- 
pendent univariate  random  variables.  Let  be  the 

rank  of  X^  among  (X^  , . . . ,XJ^) ) for  each  v = l,...,p. 
Then  corresponding  to  the  observation  matrix  X^,  we  have 
a rank  collection  matrix  R , where 


5n=  <<">>. 


i=  1 , . . . , N ; 
Vs l*...»p. 


Consider  now  p sets  of  scores  (a^  , 1 * i * N) , 1 * v * p 

generated  by  known  functions  cp  : (0,1)  -»  R in  either  of 

V 

the  following  ways: 

^i*  * *v  (i/(N+1>'  i*l#-..#»»  v=l,...,p 
aNi  = Ecpv(UN  i=l#...#N;  v-l,...,p 

where  Uj^is  the  ith  order  statistic  in  a random  sample  of 

size  N from  the  uniform  distribution  over  (0,1).  Consider 
now  the  simple  linear  rank  vector  corresponding  to  X^ 

(or  R^)  where 

S = rS(1>  s(P)l  c(v)  = r r(v)  l « * 

(SN  SN  ’ ' SN  .f^Ni  aN  (RNi  1 * " * P 

and  where  ( , l*i«N),  l*v*p  are  p-sets  of  known  HZZ 

NTIS  White  Section 

(regression)  constants.  0DC  Butt  Section 

UNANNOUNCED 

JUSTIFICATION  


D6TVIBUTI0N/AVAHA61UTY  CODES 
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In  the  next  section  we  establish  the  asymptotic  normality 
of  by  following  the  methods  of  Chernof f-Savage  (1958)  and 

Puri-Sen  (1971).  For  a different  approach,  see  H&jek  (1968) 
and  Puri-Sen  (1969).  In  Section  3 we  obtain  the  estimates  for 
the  order  of  normal  approximation  for  S„.  This  constitutes  the  gen- 
eralizations  of  the  results  of  Jureckov^-Puri  (1975),  Bergstrom-Puri 

(1977),  Puri-Rajaram  (1977),  and  Huskov^  (1977b) .where  the  problems 
are  treated  in  the  univariate  set-ups.  The  multivariate  extensions 
in  the  generality  of  our  paper  do  not  appear  to  exist  in  the  litera- 
ture so  far.  (For  a rather  special  multivariate  case,  the  reader 
is  referred  to  Huskovci  (1977a)). 

2.  Asymptotic  Normality  of  SN<  We  now  establish  the  asymptotic 
normality  of  SN  defined  in  (1.5).  We  make  the  following 


(2.1) 

(2.2) 


assumptions 

« I 

,(v) 

'Ni 


) «Pv(  t)  | * K|  (1-t)  } 


6 -a* 


|cp’  (t)  | « K{ t (1-t)  } 

V 


6 -3/2 


0<t<l 


max 

KiSH 


(2.3) 


S 

v 

where 

2 ( v) 

» = Var  S „ 

v N 


= 0 (N  j , 1 < v * P 


1 * v * P 


si1}  S<P> 


For  convenience  we  shall  take  S„  = ( ^ , where 

~N  VS,  »_  / 


the  8V  s are  given  by  (2.3). 

We  then  have  the  following  theorem. 

(v) 

Theorem  2 . 1 Let  the  scores  a„.  , lsv<P  be  defined  by  (1.3). 

Nl 


Then,  under  the  assumptions  (2.1)  and  (2.2),  for  every  vector 
X in  RP  , X.  ’ (SN -UjjJ/^’V^)*5  has  asymptotically  (as  N-»«) 


he  h(0,l)  distribution  where  u„  and 


v„  are  defined  by 
N 
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(2.4) 


(2.5) 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


A1' A*  = "A?  fvH(v,wdf2,w.iMp 

i=i  -® 

H(v)(X)  = N-1  r F^1  (X)!  F^’(x)  = P(X^tx),  1 * V *P 

vh"((^s"))  • 

~ vi  v 


8uv  = Cov  (snM)*  snV));  1 s v * p* 

Furthermore#  the  theorem  remains  true,  if  in  (2.7) 

are  given  by 

g =2  Cov  (X^^)  (X^ ) ) 

UV  i=1c'  lANi  lXNi  }'  ANi  XNi  ' U 

where 

<*>  -s  iwCw  •;«-{;! ««» 

dF^’ 


I(xsy)  = 1 if  xty,  and  0 otherwise 


Proof: 


Let  (a.(a  , ...,a  ) be  a set  of  fixed  but  arbitrary  constants. 
1 2 p 

By  the  Cramer-Wold  criterion,  it  suffices  to  show  that 


(2.10)  £ a 

v=l 


P . bN 


V N 


= U , say  is  asymptotically  normal. 


We  introduce  the  following  representation  for  Sv  . 

N 


(2.11)  H*V)  (x)  = N_1  I KxJ^  s x) 
N i_1  Nt 


5- 


H 


(v) 


(x)  is  already  defined.  (See  (2.4)  above.) 


(2.12) 

(v)  N (V) 

CN  <X>  * . E CNl 

I <X<f  S (X)) 

1 = 1 

(2.13) 

c,v) (X)  = 5 c<"> 

i-1 N1 

f^Nx)  = E[e^v)  (x)] 

We  shall  adhere  to  the  convention  of  denoting  stochastic  variables 


(hAV^  , ciV))  with  subscripts  and  non-random  functions  (H^,  C^) 
N N 


without  subscripts,  although  depending  upon  N. 


Then  the  following  inequalities  are  immediate: 
,(v)  i ^ „ max  i„(v)  i 


(2.14)  |C^  (x)  | s N |C^.;  | (x) 


(2.15)  |C(V)  (x)|  s . “ |c<">|  HWM,  1 * V p 


lsisn  1 Ni 
We  shall  use  the  representation. 


(7  16 \ c^>  - r"  m f—  u^)  d r(v)  - u(v)  + fi(v)  + B(v)  + I n(V) 

(2.16)  SN  -J  ^ Hn  (x)jdCN  (x)  - dN  + * + P2N  + .SDjN 

3=1 


( v) 

where  u is  given  by  (2.3). 
N 


(2, 

.17) 

ft(v) 

P1N 

i CD 

= 1 

J 

— CO 

<H(V> 

(X))  a 

(C(v) 
' N 

(X) 

- C(V) (x)) 

(2. 

.18) 

S(V) 

°2N 

p 00 

1 

o 

— 00 

<hhV) (x) 

-hm 

(X)) 

1 

cp 

(H(v) (x))  dc(v) (x) 

(2. 

.19) 

D<V> 

IN 

_ -1_ 
N+l 

• jv* 
— 00 

(x)  ^ 

(H(V) 

(x))  dc^v)(x) 

(2. 

.20) 

D(V) 

2N 

-r 
— 00 

<h‘V> (X) 

- H(v) 

(x)) 

cp 

V 

(H(V>(x))  d(C<V)(x)  - 

N 

C(v) 

(X)) 

(2. 

.21) 

D(V) 

3N 

j1  00 

V 

— 00 

K(sk 

HsV,<x))‘  *v<H 

<v,'x,'-(^r^v,<x' 

-H(V> 

(X)) 

dC^V) (X) 


WLm 
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Substituting  (2.16)  in  (2.10),  we  have 

(2'22)  ‘ Jl  ^ {»NV)  + (»»  + •»)  + D*  } 

The  proof  will  be  accomplished  if  we  show  the  followings 

(A)  ^ is  finite. 

/ (v)  . (vK 

P a..  \ 6 IN  ^2N  ) 

(B)  e v is  asymptotically  normal. 

v=l  8V 

p <X  3 / \ 

(C)  £ -*■  Z D^'  .»  0 in  probability. 

v=l  S j=l  IN 


Proof  of  (A) : Observe  that 

= lj  ® <H<V>  Cx) ) <JC(V,(X)| 

—CD 

s N max  |C^Y)|  f cp  (H(v)(x))  dH(v)  (x)  < 
lsi<N  W1  v 

by  assumptions  (2.1)  and  (2.15). 


+co  for  each  N 


Proof  of  (B) : For  a fixed  v , we  shall  verify  the  Liapunov  crite 


rion  for  the  normality  of  and  then  do  the  same 

P1N  P2N  v 

p 

for  the  sum  £ a (9^  + ^ an  extension  of  the  Cr“ 

\f=1  V 

inequality. 

In  fact  we  do  so  separately  for  g-^  and  9 ^ g~^  and  use  the 

IN  v P2N  v 

^-inequality. 


Since  we  are  considering  a fixed  v , we  can  drop  the  indexing 


variable  v and  simplify  our  notation  substantially 


Consider  (2.18) 

co 

$2n  = J (»N(x)  -H(x))  cp'(H(x))  dC(x) 

—«o 

Integration  by  parts  yields. 


(2.23) 

82N 

= B(x) 

|hn(x)  -H(x)j|+"  - J"  B(x)  d(HN(x)  - H (x)  ) 

—00  "“00 

where 

X 

ft 

(2.24) 

B (x) 

<p'  (H(x))  dC(x)  , xQ  arbitrary  but  H(xQ)  > 0 

xo 

We  shall  later  prove  that  the  first  term  in  (2.23)  is 

0(g)  , 8 = 2 . 

P v 


The  second  term  of  the  same  can  be  written  as 

CO 

(2.25)  j*  B (x)  d(HN(x)  -H(x))  N_1  £ -E|(y} 

—00  1—1 

We  shall  show  that 

N 

(2.26) (8N)^+e^  £ E|  8 (X^)  “ E 8 (X^)  | 2+  G -♦  0 , as  N-»oo  for  some  e > 0 . 


In  fact,  it  suffices  to  verify, 

N 

(2.27)  (8NT(2+e)  £ E|p(X  )|2+e  -»  0 

i=l 

which  by  the  and  Jensen  inequalities  implies  (2.26). 

Next,  choose  an  c>0  such  that  (2  + e)(8-Jj)>  -1  for  the 
6 given  by  the  constraints  on  . Then, 


J J 


* ,snT(2+£)  5 »a+e 


2+e 


X N2+e  max  | C | 2+6  f | f cp  ’ (H(y) ) dH(y)  | dF  (x) 
1=1  IsisN  1 ' t N1 


0 


2+e 


12 +e)  jj  ® _ 

> t max  1 11  la.,!}  f {|<p(H(x)(  | + |(p(H(x  ) )||  dF  (*) 

IsisN  1-1 


r*A 


2+e 


- 0(N  / ) f ||cp(H(x))(  + | cp  (H  (xQ)  ) | j-  dH(x)  -»  0 , as  N -♦  ® . 


because  (2+e)  (6  - Jj)  > -1  and  cp(H(xQ))  is  a constant.  Observe 

I C I 2+e 

I Nil  -1-F /2 

that  we  have  used  the  fact  max  ~ " = 0 (N  ' ) 

■%  • _ _ 2 • ® 

UuN  g 

It  remains  to  be  shown  that  in  (2.23), 


■ TK(x)  - H(x)}'  ■ °p(1> 


l6-^1!  « yN|HH(x)  -H(*)|-^|J  q) ' (H(y) ) <B(y)| 


x ,/H|iyx>  -H(x)|0<l>|  f tp' (H(y))dH  Iy)| 


K ^NlVx)  -«(x)|  |h(x)  (1  -H(x))j 


b-h 


J 


C^-inequality  yields  the  Liapunov  condition  for 
Consider  next  the  normalized  sum, 

Jl  “v  V61N  + P2N  1 • 

We  need  to  verify  the  Liapunov  condition  for  the  above 
expression. 


Let  a2  = Var(  E a > > 

\v=l  v v IN  2N 


We  shall  assume  that  a is  bounded  away  from  zero  for  all  N 
(If  not,  the  sum  is  trivially  degenerate  normal  as  N -» +»  . ) 


Write  Riw  + B 


(v)  + (v).  ^ (v) 

IN  B2N  it1  PNi 
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where  3 are  independent  random  variables  (as  already  done 
earlier) . We  have  to  show 


— 42+e)  N P 

(2.29)  lit  a .r  B e,  a “ | = 0 . 

N-»  1=1  v=l  v *v 

We  have  already  shown  that 

/•v.  * N . . 2+e 

(2.30)  lim  (§J  ; E E I eNi  I =0 

N -»  ® i=l 

(2.29)  follows  from  (2.30)  upon  noting  that  there  is  a constant 

C(p,e)  , depending  only  upon  p and  e such  that 

(v)  2+e 

•I  I “vf1  SC(P'£)  A la/6v|2+e  E'4i,|2+E 


(Generalized  C^-inequality . ) 

This  establishes  (B) . We  note  that  e depends  only  upon  6 
and  hence  the  same  choice  of  e works  for  all  v . 


Proof  of  (C) : Recall  that  we  have  to  establish 


E a S_1  1 = o (1)  . 

v-1  v v 3-1  p 


Clearly,  it  suffices  to  prove  for  a particular  v since  we  have 
a finite  sum  (1  , ...  , p)  . Again,  we  shall  drop  the  index  v • 

CO 

Consider,  I^D^l  = |- (N+iTV  1 1 J*  (x)  cp1  (H(x) ) dCN(x)  1 

S(a.T1EU'(H(xH1))1|cHi|  e v 

1=1 


i f! 
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where 


-1 


vHi  - WwtJ^n  * cK.|  . 


To  prove 
N 


Ni 


Ni 


Z V ^ -♦  0 in  probability, 
i=l 

it  suffices  to  show  that 
N 

Z Na  E|v  . Ja<  +®  , for  all  N , for  some  a , 0<a<l 

i=l 

O 

(By  particular  case  1 , page  241,  Lo6ve  (1963)). 

Take  a = 2/3  : 


8'2/3j1  eIV!2/3  iK,r2/3i!1 1 |2/3E|H(XN.)(1-H(XN.»»|  (6-3/2> 


s N1/3  max  -1 


(*1|CNi,)2/3J1  I {**<*>  <1-H(x)>} 


24. 


IsifiN 


- 1 


dF  . (x) 
Ni 


= K • 0 (1)  j*  jH(x)  (1  -H(x))|  3 dH(x)  < +«  y 


where  we  have  used  the  fact  that 


max  |5XC  |2/3=0(n~1/3)  . 
UixN  81 


Consider  next. 


00 

D2n  - ['  (HN(x)-H(x))cpMH(x))  d(CN(x)  -C(x))  . 

—CD 

By  Puri-Sen  (1971),  given  any  e>0  , 0 < 6 • < >5  , there  is  a 

constant  C(e,6')  , independent  of  N such  that  with 


probability  > 1 - c # 


13 


(2.31)  N(1+a)  E EIV  |1+a-»  0 
i=l  1 


Choose  a > 0 such  that  , (l+a)(6*-l)  > -1  ; that  is. 


0<a  < 


1-6* 


Then, 


= (l+a)  N 1-Hi 

N E E V . 

i=l  1 Ni1 

VN  max  |C  | 
* -l(l+g)|  kkN 

8 


i j,  'K'f'-Hiy}' 


(l+a)(6-i) 


-a  r"  r , (1+a)  (6*-l)  -a 

=0(l)Na  [ |H(x)  (1-H(x))J  d H(x)  = O(N^) 


which  proves  (2.31). 


Next  consider 


, • NH  (x)  N H (x) 

SD3N=S  I H N + 1 ) ~ cp(H(x) ) " £T7"  ’H(x))*'(H(x))}dCN(x) 

“CO 

We  have  to  show  that  g^D^o  in  probability. 

The  following  substitution  will  simplify  the  proof.  Let 


(2.32)  C. 


a VN 


Then  it  suffices  to  prove  C3N  * 0^(tF”) 

^Ni  —1 

Observe  that,  since  max  I I « 0(N  ) 

lii<N  ® 


we  have 


14 


(2.33)  |C3N| 


N H„(x) 


4 0(1)J  W(^  Vx)^  “ <*>(H(x))  ' (ifr  " H(x)V  (H(x))JdHN(x) 

— to 

The  proof  that  the  right  hand  side  of  (2.36)  is  o^(N  ) can 
be  found  in  Puri-Sen  (1971),  pages  401-405. 

This  proves  Theorem  2.1. 

For  a different  approach,  as  well  as  the  derivation  of  some 
tests  for  linear  hypothesis  based  on  SN  , see  Puri-Sen  (1969). 


Computation  of  V : it  is  clear  that  the  variance 

N 

contribution  is  from  the  term  + 8.!^  . It  is  easy 

* IN  w 2N 


to  check  that 


(v)  (v) 

(2.34)  g1N  + p2N 


= j1{AMV,(34Vi,)-EA'w(34i))} 


where 


ah1V)(i°  ■ ” x (c 


N |1(CNi>'  dFHj,<y)  > x0  arbitrary. 


Centering  a’  fv^(x)  at  its  expectations,  we  note  that 
Ni 


(2.35)  ,£>  * »&> 


where 


l &<&') 


*£’<*>  - h (CKi’  - c«’>  -F^l(y)}cp;<H<v)(y»dF^)(y) 

— CO 

which  yields  the  desired  approximation  (2.7). 

We  next  give  a corollary  which  extends  these  results  to  the 

case  when  the  scores  are  generated  according  to  (1.7) (b)  . 

Corollary  _2._2:  Let  m be  the  inverse  of  a distribution 

v 

* 

function  . Let  j§N  be  the  linear  rank  vector  defined  as.  in  (1.8) 

but  with  the  scores  generated  according  to  (1.4)  . Then,  under 

the  conditions  of  Theorem  2.1,  is.  a svmpto tica  1 lv  multinormal 

with  the  same  parameters  u and  V . 

N N 


multinormal 


Proof : Let 


(2.36)  S*  -I  — 

u 


* (P) 


Clearly,  it  suffices  to  orove 

(2.37)  S*(v)  - S <v)  = o (»  ) , liVip 

N N p v 

* 

which  then  entails  Sv,  - S„  -»  0 in  probability. 

~N  ~N  ~ 

Again,  since  v is  fixed,  1 s v sc  p , we  shall  drop  this  index. 
We  define 


[Nt] 

(2.38)  <pN(t)  - (gnU)  - GnU-D)  . Gn(0)  =0  , 0<t<l  , 
where  is  the  greatest  integer  not  exceeding  a . 


It  is  easy  to  check  that 
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(2.39) 


dCB(x) 


The  proof  of  Corollary  2.2  will  be  accomplished  if  we  prove 
the  following: 

Lemma  2.3:  Under  the  hypotheses  of  Theorem  2 . 1 and 

Corollary  2.2,  we  have 

(i)  lim  cp  (t)  = cp  (t) 

N-»+®  N 

(n)  )J  {<pn(^£[  hn(x))  " cp(n+T  nN(x))}d cn(xM  = °p(8) 


Proof  of  Lemma  2.3:  We  shall  only  prove  (ii)  since  the 
proof  of  (i)  can  be  found  in  pages  408-409  of  Puri-Sen  (1971) . 
By  inequality  (2.14),  we  have 


(2.40) 


“ V*>  I 


I 
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B ■ 


1 1 


But,  lira  A = 0 , by  Puri-Sen  (1971),  pages  409-411. 

H-»+»  0 0* 

S - S 
N N 

This  immediately  entails  — “ -♦  0 in  probability,  which 

O 


proves  Corollary  2.2. 

A consequence  of  Corollary  2.2  is  that  in  many  cases  of 
practical  interest  (such  as  the  normal  scores) , asymptotic 


normality  holds  with  the  same  centering  sequence  |iN  , whether 


the  scores  are  given  by  (1.3)  or  (1.4)  . Thus  theorem  2.1 

/ 

and  corollary  2.2  serve  to  unify  the  results  of  Hajek  (1968) 
and  Hoeffding  (1973),  and  also  to  some  extent  simplify  the 
results  of  the  latter  paper. 


3.  Remainder  Terms  of  S„  : In  this  section,  we  obtain 


an  estimate  for  the  remainder  in  the  normal  approximation  to  S 


Recall  that  we  can  write 


N 


(3.1)  S -u  = £ A (X  ) + D 

~N  ~N  ~Ni  Ni  ~N 


where 

(3.2) 


A (X  ) = ^a  a 

~Ni  Ni  \ Ni  Ni  ' * * *'  Ni  (XNi  ) 


(P)  ,V(P)' 


D = / n ^ ^ D 

£n  Idn  dn  ; 


(p>V 


(v)  (v)  (v)  (v) 

D v = D + D ' + D ' 
N IN  2N  3N 


livsp 


1 


Normalization  by  & is  of  no  consequence  and  was  introduced 


there  merely  to  simplify  the  details  of  the  proof.  We  shall  not 


. 


4 


use  it  here. 

We  shall  simplify  the  notation  somewhat  and  also  express 
it  in  a form  so  that  we  can  use  Corollary  17.2,  page  165  of 
Bhattacharya-Ranga  Rao  (1976) . 

We  can  write. 


<3‘3)  ~Ni  = N ANi(XNi)  ' Vi  " COV^Ni) 


Then, 


0.4)  Sb-h'*e  tn.  +dn 


-1  N 

(3.5)  V =N  e V. 

i=l  1 


Comparing  (3.5)  to  (2.7)  and  (2.8),  it  is  clear  that  the 
elements  of  V are  8 

UV 

Let,  for  x = (x^  , . . . , x^  ) , 

(3  6)  P (x)  = p/  e {1)  _ 11  (1)  * .(D  s(p)-  u(p)  < x(p)>l 

(3.6)  Fn(x)  - P(SN  uN  * x SN  UN  S X j * 

We  then  have  the  following  theorem: 

Theorem : Let  the  conditions  of  Theorem  2.1  be  satisfied. 
Let  in  addition. 


(3.7)  Sup  |cp’(t)|  = iicp'll<+®  » 1 s v s p , 

0*t<l  v v 

and  V be  positive  definite. 
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Then,  there  exists  a constant  C (p)  , independent  of  N , 
such  that. 

(3.8)  Sup  IP  (X)  - » <x>  | s C(p)  x'3/2p3 

x£1RP 

where 

(3.9)  $ (x)  is.  the  distribution  function  (p-variate)  of 

0,v  ~ 

a Gaussian  random  vector  with  mean  0 and  dispersion  matrix  V ? 

(3.10)  X = smallest  eigenvalue  of  V ; 

1 ® 3 

(3.11)  p = ” E E II  T . I!  , the  norm  being  the  euclidean  norm? 

3 N j_— Nl" 

and  -»  0 pointwise. 

Furthermore 

<3.12)  = °p(  n»vii 

Remark . Since  -*  0 , it  would  be  of  interest  to  have  an  estimate 
for  the  fluctuations  of  ^ . But  in  view  of  the  generality  of  the 
situation  and  the  fact  that  we  are  dealing  with  the  multivariate 
set  up,  it  is  not  easy  to  secure  a purely  numerical  estimate. 

(3.12)  gives  a probabilistic  bound.  The  problem  of  obtaining  a 

sharp  bound  remains  open. 

Proof:  Writing  T__=N^  E T__.  . we  have 

~N  ~Ni 

(3.13)  " Hn  = £n  + 2n  * 

Let  G be  the  distribution  function  (p-variate)  of  T . Then 

N 

(3.14)  F (X)  = G (x  - D ) . 

N ~ N ~ ~N 


Hence 
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(3.15)  |yx>  - *0  v(x)|  X |GNte-DN)  - »0<v(x  - VI 

+ l*0,v(2  ' V ■ ‘q.v1-’!  • 

We  estimate  the  two  quantities  on  the  right  hand  side 
separately. 

First,  an  application  of  Corollary  17.2  (page  172)  of 
Bhattacharya-Ranga  Rao  (1976)  yields, 

(3.16)  |Gn(x  - Dh)  - ,0_v(a  - dn>  x C(p)  ).  “3/2  p3  iT* 

uniformly  in  x , wherein  all  the  quantities  have  already  been 
defined.  ||cp']|<+«o  ensures  the  existence  of  p3 

Estimating  the  second  term  is  rather  more  involved. 

However,  since  ^ is  continuous,  the  convergence  of  the 
distribution  functions  is  uniform  and  hence 

(3'17)  S^o.v'i-V  - Vv(~*M  * an 

converges  to  zero  pointwise.  We  only  need  to  obtain  (3.12). 

Note  that,  by  the  Mean  Value  Theorem  (in  B^)  , there  is 

a point  £ on  the  Line  segment  joining  x - pN  and  x such  that 

(3.16)  *o,v(*'  -£hl<S) 

9 /*w 

where  L(t)  is  the  differential  of  at  c . 

0,v  1 

Next,  for  the  right  hand  side  of  (3.18),  we  have 


(3.19)  D^L(*) 


(v)  ± 
v=l  N 


t-  D”  (v)  #0,v(~) ' 


fix 
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Further,  for  each  v , l«v*P  , it  is  easy  to  check  that 


is  dominated  by  the  v marginal  density  . 


Hence , | 


. (v)  1 V 


Consequently, 

(3.19)  |4H|  * Z »'1|DnV>I  (2'’)"'5  • 

V=1  v 

Thus,  (3.12)  will  be  established  if  we  show  that  for  each  v , 

(u)  . N (v) 

(3.20)  V V'  . 0p(||¥v||  I |CN.|)  . 

Again,  we  shall  drop  the  index  v in  the  rest  of  the  proof. 

Rearranging  terms  of  D„  , we  get, 

N 

(3.2D  dh  = J*{<p(575;  V*))  ' 

—CD 

CD 

- I*  («N(x)  - H(x))cp'  (H(x))  dC(x)  . 


In  order  to  simplify  the  proof,  we  shall  drop  the  factor 
— . In  view  of  the  fact  that  cp1  is  bounded,  the  conclusion 
will  not  be  affected. 

Consider  <p(HN(x))  - cp(H(x))  . By  the  Mean  Value  Theorem, 
Cttyx))  - cp (H(x) ) = (Hjj(x)  - H(x))  cp * (§jj(x)) 

for  some  J <*>  • 
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Hence, 


(3.22)  |„/  N j (Hn(x)  - H(x))cp  ($N(x) ) d CN(x)  \ 

“W 

00 

* fc'll  J y»|H„(x)  - H(x)  |d|CN|  <x) 


Let  e >0  be  given.  Then  by  Puri-Sen  (1971),  there  is  a 
constant  C(e)  independent  of  N such  that  with  probability  >1-$ 


J N | Hjj(x)  - H (x)  | < C(e) 


Thus  in  (3.22),  with  probability  > 1-e  , 


(3.23)  |pH(HN(x)-H(x))?'(|N(x))dCN(x)|  * C (e)  \\cp*  \\  ^1^1 


The  proof  of 


(3.24)  \\  JV  (Hn(x)  -H(x))<p'  (H(x))dC(x)|  £ C(g)||«p,||  S \C^\ 

— 00 

in  probability  is  identical. 

(3.23)  and  (3.24)  establish  (3.20)  and  (3.19),  which  in  turn 
entail  (3.12).  The  proof  is  completed. 


Remark:  From  the  expressions  for  S and  g , it  is 
(JV  v 

clear  that  \ is  extremely  hard  to  compute.  Recall  that  X is 
the  smallest  eigenvalue  of  the  dispersion  matrix  V . However, 
if  the  components  of  the  vectors  X^  satisfy  the 
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condition  of  "weak  dependence"  in  the  sense 

(3.25)  a2  - s Is  ,|>0  , v-l,...,p  ; 

v U*u' 

the  classical  Grischgorin  theorem  (cf.  Dahlquist-B jorek,  (1974)), 
permits  us  to  replace  x by  much  simpler  expression  (3.25)  , 
assuming  without  loss  of  generality,  x>0  • if  one  is  prepared 
to  work  numerically,  much  better  estimates  for  X can  be  obtained 
(op.cit.  Dahlquist-B jorck) . 


I 
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